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ABSTRACT 

A  two-dimensional,  vertical-plane,  computer  simulation  model  de¬ 
scribing  the  dynamics  of  a  wheeled  vehicle  exiting  from  the  riverine 
environment  is  presented.  The  model  includes  the  effects  of  soft  soil, 
suspensions,  buoyancy,  drive  train  characteristics  and  the  inertial 
reactions  of  the  unsprung  masses.  The  body  is  modeled  as  a  composite  of 
rectangular  sections;  the  tires  are  rigid.  The  stream  bottom  and  bank 
surface  can  be  any  specified  arbitrary  geometry;  their  soil  properties 
are  specified  by  the  land  locomotion  soil  value  system. 
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I.  INTRODUCTION 

The  philosophy  and  initial  efforts  in  simulating  the  exiting  of  an 

off-road  vehicle  from  a  body  of  water  was  presented  by  Sloss,  Ehrlich  and 

1  2  3 
Worden  in  the  second  volume  of  a  series  of  reports  ’  on  general  studies 

of  vehicles  in  the  riverine  environment.  These  reports  defined  the  total 

river-crossing  problem  as  being  composed  of  ingress,  water  performance  and 

egress.  The  egress  portion,  being  the  most  difficult,  was  considered  in 

need  of  the  most  urgent  analytic  study. 

In  response  to  this  need,  a  simple,  vertical  plane  computer  simu¬ 
lation  model  was  written  and  validated  with  a  small  vehicle  in  the  idealized 
environment  of  a  laboratory  river  simulator.  This  model  was  restricted  to 
rigid,  planar  banks,  consisting  of  only  a  single  ramp.  The  surface-wheel 
interaction  was  modeled  by  a  single  spring-damper  wheel  deflection  mechan¬ 
ism,  the  tire-surface  coefficient  of  friction  and  the  bank  slope. 

4 

Muddappa  and  Baker  attempted  to  extend  this  model  by  including  more 
general  soil-wheel  interactions  based  on  the  land  locomotion  soil  value 
system.  However,  since  this  model  employed  only  equilibrium  equations, 
it  could  not  truly  simulate  dynamic  situations. 

The  model  described  herein  is,  basically,  an  extension  of  the  Sloss- 
Ehrl ich-Worden  model.  It  extends  the  simulation  to  arbitrarily  shaped 
banks  with  arbitrary  soil  parameters,  but  it  restricts  the  wheel  to  being 
solid.  The  inclusion  of  flexible  wheels  awaits  the  development  of  flexible 
wheel  models  which  are  capable  of  simulating  travel  over  both  hard,  rough 
terrain  and  soft  soils. 

Since  it  is  anticipated  that  this  simulation  will  be  used  to  describe 
motion  over  both  obstacles  and  hard  terrain,  a  much  more  complete  set  of 
motion  equations,  including  many  inertial  terms  ignored  in  the  earlier 
models,  are  included. 
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A  major  feature  required  for  simulation  of  wheel-soil  interaction 
is  the  calculation  of  wheel  slip.  This  required  the  complete  dynamic 
modeling  of  the  forces  tending  to  spin  the  wheels,  including  those  from 
the  engine-drive  train,  soil  shear  resistance,  and  hydrodynamic  drag  on 
the  spinning  submerged  tire.  Thus  the  model  presented  here  represents  a 
major  increase  in  sophistication  over  previous  efforts. 
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II.  MODEL  DESCRIPTION 


The  Vehicle 

The  vehicle  is  treated  as  a  multi-mass  system  consisting  of  a  body 
and  any  number  of  axles/'  The  entire  vehicle  is  treated  symmetrically 
in  the  lateral  direction  and  is  therefore  assumed  restricted  in  roll, 
sway  and  yaw.  This  results  in  a  vertical  plane  simulation  allowing  only 
three  degrees  of  freedom:  surge  (x) ,  heave  (z)  and  pitch  (9).  Figure  1 
presents  the  coordinate  system  employed,  using  a  two-axle  vehicle  as  an 
example. 

The  axles  are  assumed  to  be  thin  rods  with  their  associated  unsprung 
mass  concentrated  on  a  line  connecting  the  left  and  right  wheels.  The 
wheels,  attached  to  the  axles  at  each  end,  are  assumed  to  have  rectangular 
cross-sections  and  are  allowed  only  to  translate  in  the  vehicle  heave  (z) 
direction.  All  wheels  are  constrained  to  rotate  together  so  that  only 
one  wheel  rotational  degree  of  freedom  is  allowed;  thus  the  vehicle  is 
assumed  to  have  all-wheel  drive  with  no  differential  between  front  and 
rear  axles.  Due  to  the  left-right  assumption  of  symmetry,  both  wheels  of 
an  axle  rotate  at  the  same  speed. 

The  vehicle  body  is  considered  to  consist  of  a  number''  of  rectangular 
boxes  with  rectangular  cross  section.  These  boxes  and  the  wheel  assemblies 
generate  buoyant  forces  when  submerged.  For  linear  motion  the  entire  mass 
of  the  body  is  assumed  to  be  concentrated  at  its  center  of  gravity. 

The  axles  are  joined  to  the  body  by  a  spring,  damper  and  bump-stop 
system  (see  Figure  2).  Multiple  rates  for  the  spring  and  both  bump  stops 
are  allowed.  The  damper  is  assumed  to  have  different  rates  and  blow-off 
forces  for  jounce  and  rebound,  but  to  be  only  velocity  dependent. 

-'Currently,  for  convenience,  the  simulation  is  restricted  to  8  axles  or  a 
1 6x  1 6  vehicle;  however,  this  restriction  is  easily  removed  without  com¬ 
promising  the  integrity  of  the  model. 

Currently,  up  to  ten. 
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The  drive  train  dynamics  are  governed  by  a  curve  relating  engine 
RPM  to 'torque  at  full  throttle.  The  transmission  is  allowed  up  to  nine 
different  gear  ratios.  A  single  overall  final  drive  ratio  is  assumed. 

The  Driver 

The  "driver  model"  currently  assumes  full-throttle  at  the  lowest 
gear  ratio  possible  for  the  required  torque  within  the  gear  range  speci¬ 
fied  at  the  commencement  of  each  simulation  run.  The  best  driver  strategy 
for  egress  is  not  yet  known;  hence  it  is  expected  that  this  portion  of  the 
model  may  eventually  undergo  many  substantial  revisions. 

The  Bank 

The  bank  is  described  by  a  table  relating  bank  surface  elevation 
and  distance  relative  to  the  water/bank  intersection  (see  Figure  3). 
Elevation  and  soil  parameters  may  change  anywhere  along  the  bank. 

The  bank  coordinates  are  used  as  the  earth-fixed  frame  of  reference; 
vehicle  position  is  specified  in  bank  coordinates.  The  origin  of  this 
reference  frame  is  at  the  water  bank  intersection;  the  horizontal  water 
surface  is  the  negative  x’-axis;  the  z'-axis  is  vertical,  positive  down. 
The  bank  contour  is  specified  by  a  number  of  convenient  points,  P.  ,  (see 
Figure  3)  from  which  an  auxiliary  program,  BNKPR0.F4,  described  in 
Reference  5,  constructs  the  detailed  bank  profile  by  fitting  a  cubic 
equation  to  four  consecutive  straddling  points,  with  appropriate  modifi¬ 
cations  near  the  end  points. 

Soil  strength  characteristics  along  the  bank  are  specified  by  a 

set  of  x'-coordinates  (S .  in  Figure  3).  The  region  immediately  to  the 

left  of  each  boundary  point,  S.  ,  is  assumed  to  have  homogeneous  soil 

characteristics  as  specified  by  the  land  locomotion  soil  value  system 

parameters.  The  particular  manner  in  which  they  are  used  here  is  given 
z 

by  Hiszklevitz. 
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FIGURE  3.  TYPICAL  BANK  PROF  I 
P.  ARE  BANK  PROFILE  DATA  POIN 
ARE  BOUNbARIES  BETWEEN  HOMOGENEOUS  Si 
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III.  EQUATIONS 

The  state  variables  of  the  vehicular  equations  of  motion  are  both 
the  position  and  the  velocities  of  the  various  vehicular  components. 
Simplification  of  calculation  is  achieved  by  writing  the  differential 
equations  of  motion  in  the  vehicle-based  coordinate  system  shown  in 
Figure  1  and: 

1.  Calculating  the  unbalanced  forces  resulting  from  the 
velocities  and  positions  of  each  vehicular  component, 

2.  Converting  the  unbalanced  forces  of  Step  1  into  resultant 
acce lerat ions , 

3.  Integrating  the  appropriate  accelerations  of  Step  2  and 
velocities  of  Step  1  to  achieve  new  velocity  and  positional 
data, 

4.  Converting  the  new  vehicular  coordinate  positions  and 
velocities  of  Step  3  to  earth-fixed  coordinates, 

5.  Repeating  the  cycle  by  returning  to  Step  1  using  the 
velocities  and  positions  calculated  in  Step  4. 

The  fact  that  individual  accelerations  and  forces  are  not  calculated  from 
a  balance  of  forces  means  that  no  equilibrium  analysis  is  employed. 

For  example,  traction  at  each  wheel  is  not  computed  by  calculating 
the  net  DBP  required  to  overcome  various  resistances  and  allocating  DBP 
to  the  various  wheels  by  some  method.  Rather,  individual  wheel  unbalances 
and  interferences  with,  say,  the  soil  are  used  to  calculate  forces  and 
accelerations  which,  when  converted  to  actual  motions,  will  redress  the 
unbalances.  This  general  approach  will  be  given  explicit  form  in  this 
section. 

Forces  governing  vehicular  and  component  motion  result  from  the 
following  convenient,  though  not  mutually  exclusive,  categorization: 
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o  soil/wheel  compaction  and  shear 
o  water  buoyancy 

o  viscous  drag 

o  drive  train  thrust 

o  springs  and  dampers 

o  gravity 

Actually,  the  only  two  sources  of  forces  are  gravity  and  the  energy  stored 
in  the  fuel  which  generates  the  drive  train  thrust.  However,  for  con¬ 
venience,  the  results  of  their  application  may  be  categorized  as  described 
above,  if  care  is  taken  not  to  include  a  particular  force  under  more  than 
one  category.  The  equations  describing  these  forces  in  the  vehicular 
coordinate  system  are  presented  below. 

Soil -Wheel  Interactions 

The  soil -wheel  model  employed  here  has  been  reported  in  detail  by 
Mi  szklevi tz.^  Only  a  brief  outline  of  this  approach,  along  with  the 
resulting  equations,  will  be  given  here. 

The  major  forces  resulting  from  the  soil-wheel  interactions  are  the 
lift  (or  support)  forces,  primarily  as  a  result  of  soil  compaction  due  to 
sinkage  resistance  and  the  driving  forces,  primarily  the  result  of  soil 
shear  due  to  the  rotation  of  the  wheel  caused  by  drive  train  torque.  In 
addition,  the  driving  shear  force  also  has  a  lift  component  and  the  com¬ 
paction  gives  rise  to  a  drag. 

Due  to  the  simplification  of  the  vehicle  suspension  (motion  of  the 
wheel  in  the  z-direction  only)  the  pitch  moments  about  the  vehicle  CG 
generated  by  the  four  forces  mentioned  above  are  treated  somewhat  different 
ly  depending  on  the  coordinate  direction  in  which  they  occur.  Forces  on 
the  wheel  in  the  z-direction  result  in  motion  of  the  wheel  relative  to  the 
body.  This  relative  motion  results  in  suspension  forces  which  are  applied 
to  the  body  at  suspension  attachment  locations.  Thus  no  pitch  torques  or 
heave  forces  resulting  from  soil-wheel  z-direction  forces  are  applied  to 
the  body  directly;  they  will  eventually  appear  as  suspension  forces  and 
moments  on  the  body.  However,  no  relative  x-direction  motion  between  the 


9 


wheel  and  the  body  is  allowed;  therefore,  soil-wheel  x-direction  forces 
give  rise  to  pitch  torques  and  surge  forces  on  the  body  directly.  The 
moment  arms  of  these  so i 1  - insp i red  pitch  torques  are  the  z-di stances  from 
the  CG  of  the  body  to  the  location  in  the  soil -wheel  contact  patch  where 
the  soil  forces  are  generated  (see  Figure  4). 


FIGURE  4.  SCHEMATIC  SHOWING 
SOIL-GENERATED  FORCES  AND  MOMENTS. 
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At  any  instant  in  time,  the  wheel  under  consideration  has  a  position 
relative  to  the  bank  which  may  or  may  not  cause  soil-wheel  interference. 

If  there  is  no  interference,  there  are  no  soil-wheel  forces.  Interference 
is  detected  by  calculating  the  z'-coordinate  of  many  locations  on  the 
wheel  circumference  and  comparing  that  with  the  corresponding  z 1 -coord inates 
of  the  bank.  When  interference  is  detected,  the  locations  of  entry  and 
exit  interference  with  the  original  and  current  bank  contours  are  deter¬ 
mined  and  expressed  as  angles  of  the  wheel  relative  to  the  horizontal 
(see  Figure  5)* 


FIGURE  5-  SCHEMATIC  DEFINING  SOIL-WHEEL 
ENTRY  AND  EXIT  LOCATIONS. 
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The  actual  ground  contact  patch  extends  only  from  the  entry  point, 

Y|  ,  to  the  exit  point,  Y2  ’  however,  the  bank  originally  had  a  contour 

different  from  the  current  position,  the  soil  has  experienced  some  pre¬ 
compaction.  Therefore,  to  calculate  the  soil  forces  on  the  wheel,  the 
normal  stress  distribution  is  computed  as  if  soil-wheel  contact  extended 
from  Y]0  to  Y2  >  but  is  integrated  only  from  Yj  to  Y2  to  calculate  total 
normal  force.  (See  Figure  6.)  Support  for  this  method  of  calculating 
lift  forces  is  found  in  Bekker^  and  elsewhere;  the  entire  argument  is  pre- 
sented  in  detail  by  Miszklevitz. 


The  normal  force  developed  in  the  contact  patch  depends  primarily 
on  the  maximum  sinkage  and  the  stress  distribution.  Here  a  symmetric 
quadratic  stress  distribution  was  postulated,  not  because  it  is  known  to 
be  true,  but  because  the  deviations  from  measured  distributions  are  not 
great  and  the  resulting  normal  stress  equation  can  be  analytically 
integrated.  If  a  is  the  maximum  normal  stress  and  Y  is  any  angle 
between  Yjq  and  y2  >  then  the  normal  stress  at  Y  is  postulated  to  be 
given  by 


whe  re 


2 

ct(y)  =  AY  +  By  +  C  (1) 

ho 

_ max _ 

^  Y1 0“ Y2^  (Y2-Y10^ 

B  =  -(y10+y2)a 

c  =  Y10Y2A  (2) 


This  results  in  a  distribution  shown  in  Figure  7.  The  maximum  normal 

stress,  a  ,  is  calculated  from  Bekker's  equation 
’  max  ’ 


a 

max 


(3) 
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FIGURE  6.  NORMAL  STRESS  DISTRIBUTION  UNDER  A  REAR  WHEEL 
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Y 
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FIGURE  7-  ASSUMED  NORMAL  STRESS  DISTRIBUTION 


SOIL  PROFILE 


FIGURE  8.  LOCAL  GROUND  SLOPE  AND  SINKAGE  FOR  A  RIGID  FRONT  WHEEL 
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where  z^  is  the  sinkage  at  the  center  of  the  contact  patch  (not  necessarily 
the  maximum  sinkage,  see  Figure  6)  which,  for  a  wheel  of  radius  r  ,  is 
given  by: 

zm  =  r[sin(0G - -)  -  s  in ( eG-^10)  1  (*0 

Here  0  is  the  local  ground  slope,  being  the  ang le  wh ich the  horizontal  of 

U 

the  chord  connecting  the  intersections  of  the  wheel  makes  with  the  original 
bank  surface,  Yjq  and  Y20  (see  Figure  8).  Also  kc,k^  and  n  in  Equation  (3) 
are  Bekker's  soil  strength  parameters. 

The  soil  shear  stress-strain  distribution  in  the  contact  patch  was 

8 

proposed  by  Janos i  and  Hanamoto  as 

t(y)  =  (c+o(y)  tan?)  (l-e“j/k)  (5) 

Here,  c  and  cp  are  the  Coulomb  soil  shear  strength  parameters  and  j  is  the 

q 

soil  deformation,  which  is  given  (after  Schuring  and  Belsdorf  )  by 

j  =  r  [Yj-Y+(  1-s)  (cos  (Y-  Bq)  -  cos  (Y^  9g))  3  (6) 

where 

rq  -  u„ 

g  t  \ 

s  _  -  _  ^|ie  longitudinal  slip  of  the  wheel  (7) 

rqw 

q^  =  the  angular  velocity  of  the  wheel 

Ug  =  the  velocity  component  of  the  wheel  center 

in  the  9r  direction 
b 

The  above  equations  for  o(y)  and  t(y)  are  used  to  assign  normal  and 
shear  stresses  to  each  point  along  the  wheel  contact  patch  between  Yj  and 
Y2  •  These  stresses  are  converted  to  vehicle  coordinates  and  integrated 
over  the  contact  patch  to  yield: 
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But : 


Therefore: 


i 

Compression  drag  =  F  =  -  br  f  o(y)  cos(y~9)  d Y 


V, 


Compression  lift  =  F  =  br  J  o(y)  sin(Y-9)  dY 


Vo 


i 

Shear  thrust  =  =  -  br  J  t(y)  sin(Y-O)  dY 


Vo 


Y 


Shear  lift  =  F  =  -  br  J  t(y)  cos(y-9)  dY 

Z  Vo 


Total  shear  moment  on  wheel  =  T  =  br^  J  t(y)  dY 


Compression  drag  moment  about  vehicle  CG  = 

,Y1 

Mxa.  =  br  J  (z.+  r  cosa)  o(y)  cos(y-B)  dY 


Shear  thrust  moment  about  vehicle  CG  = 


I 

Mxt  =  br  J  (z.+r  cosa)  t(y)  sin(Y-9)  dY 


TT 

ot  =  2  +  Y-9 


Y 


M  =  z.F  +  br  r  ct(y)  sin(Y-S)  cos(y-Q)  d 

xo  i  xo  J  v  v  '  v  ' 

\ 


(8) 

(9) 

(10) 


(11) 


(12) 


(13) 


(14) 


(15) 
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Y 


Mxt  =  ziFxT  +  br  J  T(Y)  sin(Y-e)  cos(Y-e)  dY 

Y2 


(16) 


The  total  force  parallel  to  the  vehicle  x-axis  is  given  by: 


F  =  F  +  F 
sx  xt  xa 


(17) 


and  that  parallel  to  the  vehicle  z-ax?s  is: 


F 

sz 


+  F 


zo 


Also: 


M  =  M  + 
sx  xo 


XT 


(18) 


(19) 


The  equations  for  F  ,  F  ,  and  M  were  integrated  analytically  by 
^  g  xcr  zcr  xa 

Miszklev i tz.  The  equations  for  F  ,  F  and  M  are  not  readily  inte- 

XT  Z  T  XI 

gratable  analytically;  so  they  are  integrated  numerically  whenever  needed. 
Details  are  presented  in  the  Section  IV  describing  the  soil  wheel  inter¬ 
action  subroutine,  SOILF. 

Wheel  Buoyancy 

For  purposes  of  the  calculation  of  the  buoyant  forces  on  the  wheels, 
it  was  assumed  that  the  wheels  were  toroids  of  rectangular  cross-section 
mounted  on  thin  rims  which  do  not  displace  a  significant  amount  of  water. 

For  a  given  z'-position  of  the  center  of  the  wheel,  the  area  of  the 
wheel  torus  under  water  multiplied  by  the  wheel  width  would  give  the 
volume  of  water  displaced.  Multiplication  of  the  displaced  volume  by  the 
density  of  water  yields  the  buoyant  force  on  the  wheel.  Thus,  for  the 
disc  defined  in  Figure  9: 

0  =  cos  ^  (z'/r) 

w  \  '  y 
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2  2 

the  area  of  wheel  swept  by  2  9  =  TTr  (9  /tt)  =  r  9 

r  7  w  v  w  7  w 


1  [2  2  / 2 

the  area  of  triangle  underwater  =  y  z  1  2//r  -z 1  =  z1  vr  -z1 


2  2  /  2  2 

Thus  the  total  area  submerged  =  Ag  =  TTr  -  r  9  +  z1  Jr  -z1  (20) 


FIGURE  9.  SIGNIFICANT  VARIABLES  IN  THE  WHEEL  BUOYANCY  CALCULATIONS 

In  Equation  (20),  r  is  either  the  radius  of  the  wheel,  r^  ,  or  the  radius 
of  the  rim,  r^. 

Then,  the  wheel  buoyant  force  is 

F,  =  p  b(A  -  A  D)  (21) 

b  '  sw  sR7  x  7 

The  line  of  action  of  the  wheel  buoyant  force  is  always  considered  to  be 
in  the  vertical  (z 1 -negat ive)  direction  at  the  wheel  center. 

The  calculations  described  above  for  both  the  wheel  and  rim  buoyant 
forces  are  made  for  each  wheel  in  the  subroutine  UPDATE,  described  in 
Section  IV. 

Wheel  Rotational  Drag 

g 

From  a  study  by  Ehrlich  and  Nuttall,  it  is  known  that  the  hydraulic 
rotational  drag  on  a  spinning  tires  is  sufficiently  high  to  consume  com- 
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pletely  the  torque  derived  from  the  M151  engine  at  approximately  25  mph. 
Data  for  the  performance  of  an  M 1 5 1 ,  with  wheels  fully  submerged  is  given 
in  Table  1 . 

Figure  10  shows  the  estimated  relationship  between  wheel  rotational 
drag,  D  (expressed  in  ft-lb),  and  wheel  speed,  qw  (expressed  in  rpm) , 
for  the  standard  7.00x16  NDCC  tire  mounted  on  the  M 1 5 1  1/4-ton  truck.  The 
curve  fitted  to  the  data  of  Table  1  is: 

0„  =  a(qw)b  (22) 

where:  a  =  .008123 

b  =  1.834 


The  speedometer  reading,  not  the  linear  speed  of  the  vehicle. 
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If  is  expressed  in  radians  per  second,  a  =  6.112  for  Equation  (22). 

This  relationship  is  currently  in  the  simulation  in  the  subroutine  UPDATE, 
described  in  Section  IV.  A  change  to  other  size  wheels  and  tires  may 
require  a  change  to  the  program. 

Driving  Forces 

For  the  purposes  of  this  simulation,  it  is  assumed  that  all  wheels 
are  rigidly  connected  to  the  drive  train  with  no  differential  between 
adjacent  wheels  or  between  front  and  rear  axles;  the  vehicle  has  a  mechani¬ 
cal  transmission;  all  final  drive  ratios  are  assumed  equal;  and  the  throttle 
is  assumed  fully  open.  The  parameters  describing  such  a  drive  train  are: 

o  the  plot  of  engine  RPM  vs.  torque  at  wide  open  throttle. 

o  the  number  of  gears,  Nr  ,  in  the  transmission  and  the 

b 

ratio  of  each  gear,  G.  i=l,  ...»  N„  . 

I  u 

o  the  final  drive  ratio,  Gp  ,  which  includes  all  gearing 
from  the  transmission  to  the  wheels. 

o  the  efficiency  of  the  transmission  gears,  T|G  ,  and  the 
final  drive  gear  train,  T|  . 

o  the  engine  speed,  ,  at  or  below  which  the  next  lower 
gear  is  selected  by  the  driver  while  decelerating. 

o  the  engine  speed,  ,  above  which  the  next  higher  gear 
is  selected  by  the  driver  while  accelerating. 

Given  the  wheel  speed,  q^  ,  the  engine  speed  is  calculated  for  each  trans¬ 
mission  gear  setting  by: 

^Ri  =  Gi  GF  qw  ^or  '  =  Ng  (2^) 

The  lowest  value  of  i  for  which  SD  <  E  .  ^  UD  is  the  selected  gear  ratio. 

K  K I  K 

For  that  gear  the  engine  torque,  Tc  ,  at  the  engine  speed,  E  .  ,  is 

t  K  l 

calculated  by  interpolation  from  the  table  relating  engine  speed  to  engine 
torque  at  wide  open  throttle.  The  torque  available  to  the  wheels,  T^  , 
is  then  calculated  by: 
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Tw  -  TE  Gi  T|G  GF  %  (24) 

These  calculations  are  done  in  the  subroutine  UPDATE,  described  in 
Section  IV. 

Suspension  Forces 

Figure  2  defines  the  parameters  of  the  suspension.  Since  the  simu¬ 
lation  is  restricted  to  the  vertical  plane,  the  complicated  suspension 
motions  resulting  in  lateral  side  slip,  camber  and  track  variations  are 
not  included  here.  Wheel  base  variations  were  considered  negligible. 

These  considerations  led  to  simple  wheel  motions  restricted  to  be 
only  in  heave  (parallel  to  the  vehicle  z-di rection) .  A  no-load  spring 
extension  distance,  z^  ,  is  postulated  which  allows  spring  forces,  , 
to  be  calculated  by 

S  =  (z  -z)  k  for  z  ^  z  s  z 
p  R  '  u  a 

Sp  -  <2lC2u>  k  +  <Vz)  ku  for  2  <  2u 

sp  “  (ziC2*>  k  •  <2-2.e>  ki  for  zi  <  2  <25) 

The  convention  used  here  is  that  the  suspension  forces  which  act  to 
separate  the  wheel  from  the  body  are  positive. 

Damper  forces  depend  on  the  vertical  velocity  of  the  wheel,  z  ,  and 
always  oppose  the  motion  of  the  wheel.  Most  hydraulic  shock  absorbers 
are  constructed  to  give  lower  damping  rates  for  compression  (upward  move¬ 
ment)  than  for  rebound.  In  the  current  model  there  is  provision  for  one 

damping  rate,  d  ,  for  negative  values  of  z  and  another,  d  ,  for  positive 
c  r* 

values.  Furthermore,  hydraulic  shock  absorbers  have  "blow  off"  characteris¬ 
tics;  that  is,  they  are  capable  of  delivering  forces  only  up  to  a  given 
maximum  value,  F,  ,  usually  different  in  each  direction,  F,  and  F, 

Thus,  the  suspension  forces  due  to  hydraulic  dampers  are: 
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=  maximum 

'V-  -Fdc* 

for 

z  <  0 

=  minimum 

ldr*-  Fdr> 

for 

O 

A1 

N 

The  total  suspension  force  is  given  by 


(26) 


(27) 


The  above  calculations  are  included  in  the  subroutine  UPDATE,  described 
in  Section  IV. 

Body  Buoyant  Forces 

For  purposes  of  buoyant  force  and  center  of  buoyancy  calculations, 
the  body  is  divided  (in  the  vertical  plane)  into  rectangular  boxes  whose 
sides  are  parallel  to  the  vehicle  (x-z)  coordinate  axis.  Each  corner  of 
each  box  are  specified  by  one  x-  and  one  z-coord inate .  The  box  is  assumed 
centered  along  the  y-axis  and  is  assigned  only  a  width,  wg.  . 

The  equation  of  the  waterline  in  the  vehicle  coordinates  is  given  by 

z  =  s  x  +  z  (28) 

w  w  v  / 

where:  the  slope  of  water  line  with  respect  to  the  vehicle  is 

s  =  tan  0 
w 

the  z-intercept  of  waterline  with  the  vehicle  coordinate  is 

z  =  -z‘  /cos  9 
w  CG 

TT  TT 

The  pitch  angle  9  ,  is  restricted  to  be  between  -  j  and  y  .  The  simulation 
ends  if  0  goes  out  of  that  range,  since  the  vehicle  has  expe rienced  a  capsize. 

Required  for  the  calculations  of  the  center  of  buoyancy  are  the  equa¬ 
tions  which  locate  the  centroid  of  a  triangle,  given  the  x-  and  z-coordi- 
nates  of  each  apex.  Thus,  if  the  triangle  has  coordinates  given  in 
Figure  11,  the  centroid,  c  ,  has  coordinates 

x1+x2+x3  z1+z2+z3^ 

V  3  ’  3  J 
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FIGURE  11.  CENTROID  OF  TRIANGLE 


In  general,  the  buoyancy  is  calculated  for  each  box  by  calculating 
the  area  of  the  box  below  the  waterline,  then  multiplying  by  the  width 
of  the  box  and  the  density  of  water.  The  center  of  buoyancy  is  the 
centroid  of  the  area  under  water.  The  line  of  action  of  the  buoyant 
force  is  in  the  vertical  upward  (negative  z1)  direction  through  the 
centroid. 


Given  the  vehicle  location,  x^  and  z^  ,  and  the  attitude,  9  , 
the  z 1 -coord i nates  of  each  of  the  box  corners,  z'  ,  can  be  calculated  by 


Bj 


zBj  =  ZCG  -  XBj  5in  6+  ZBj  cos 


Here  the  subscripts  are  coded  as: 

subscri pt 

j  =  1 
j  =  2 
j  =  3 

j  =  4 


j  =  1, 2, 3, 4 


corner  of  box 

lower  front 
lower  rear 
upper  front 
upper  rear 


(29) 


The  values  of  z^  are  then  compared  to  the  water  level  (always  at  z1  =  0) 
and  the  various  possible  positions  of  the  box  are  categorized  into  the 
following  cases : 
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Case  0.  Box  entirely  out  of  water: 


=  the  force  of  buoyancy  =  0 

=  the  x-coordinate  of  the  centroid  =  undefined 
=  the  z- coordinate  of  the  centroid  =  undefined 


(30) 


Case  1.  Only  the  lower  front  corner  ( j  =  1 )  under  water  (see 
Figure  12) : 

If  the  pitch  angle,  181,  is  less  than  .01,  it  is  assumed 
that  the  entire  box  is  out  of  the  water  and  the  calculations  of  Case  0 
are  performed. 


k 


FIGURE  12.  BUOYANCY  CALCULATION,  CASE  1 


and 


For  1 0 1  s  .01,  the  coordinates  of  A  are  (x^,  z^)  where 

the  coordinates  of  B  are  (xD1,zD)  where  zD  =  x„,  s  +  s  . 

BIB  B  B  I  W  W 

The  area  submerged  is 


*A  ~ 


s 


w 


9 


A 

s 


*A}  (zBl  "  ZB); 


the  force  of  buoyancy  is 

fb  =  p  as  wB  ; 


(30 
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and  the  centroid  is  located  at 


x 


FB 


2x 


B 1 


+  x. 


z 


FB 


(31  cont'd) 


Case  2.  Only  the  lower  rear  corner  (j  =2)  is  under  water  (see 
Fi gure  13)  : 

If  101  ^  .01,  it  is  assumed  the  entire  box  is  out  of  the 
water  and  the  calculations  of  Case  0  are  performed. 


3 


FIGURE  13.  BUOYANCY  CALCULATION,  CASE  2 

For  1 9 1  2:  .01,  the  coordinates  of  A  are  (x.,  zD„)  where  x 

A  d  L  * 

and  the  coordinates  of  B  are  (xD„,  zD)  where  z_  =  xD„  s  +  z  • 

oZ  d  B  dz  w  w 

The  area  submerged  is 

As  =  2  ^ZB2  "  ZB^  (*A  “  *B2^  ’ 


'B2 


w 


w 
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the  force  of  buoyancy  is 


FB  -  P  As  WB  ; 

and  the  centroid  is  located  at 

2xb2  +  xA 
XFB  "  3 

2zB2  +  ZB 
ZFB  3 


Case  3.  Both  lower  corners  (j  =  1,  j  =  2)  only  are  under  water 
(see  Figure  14) : 

3  4 


FIGURE  14.  BUOYANCY  CALCULATION,  CASE  3 

The  coordinates  of  A  are  (xB2>  zft)  where  z A  =  xg2  sw  +  zw  ,  the 
coordinates  of  B  are  zg)  where  z^  =  Xgj  sw  +  zw  ,  and  the  coordi 

nates  of  C  are  (x_,  z£)  where 

XC  =  xB2’  zC  =  ZB  lf  ZB  ^  ZA 
XC  ~  XB 1  ’  ZC  =  ZA  ‘f  ZB  <  ZA 
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Let 


zo  ■  ZA  for  zb  *  za  - 


ZD  *  ZB  for  ZB  <  ZA 


A^  =  area  of  triangle,  ABC  , 

Ad  =  area  of  rectangle  12CB  for  zD  ^  z.  ,  or  12AC  for  zn  <  z.  , 
R  3  BA*  BA’ 

1  =  XB1  '  XB2 


Then 


A, 


T  2  ^ZC  "  ZD^ 


and  the  centroids  of  the  area  A^  are 


XT  = 


XB1  +  XC  +  XB2 


ZA  +  ZB  +  ZC 


A 1  so 


ZT  ■ 


AR  "  *<ZB2  -  ZC> 


and  the  centroids  of  the  area  AD  are 

K 


XR  = 


XB1  +  XB2 


ZR  = 


ZB2  +  ZC 


The  submerged  area  is  then 


A  —  At  +  A  ; 
s  T  r  ’ 


the  force  of  buoyancy  is 
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FB  =  PAs  WB  5 

and  the  centroid  of  the  submerged  area  is  located  at 

Vt  +  xrar 

XFB  A 


(33) 


ztat  +  zrar 


'FB 


Case  4.  Both  front  corners  (j  =  1,  j  =3)  only  are  submerged 
(see  Figure  15) : 

4 


TT 


Since  9  is  never  less  than  -  j  ,  C  can  be  assigned  without  ambiguity 
e  coordinate: 
respectively,  and: 


and  the  coordinates  of  A,  B  and  C  are  (x^,  2gj),  (XB>  z^)  and  (x^,  z^), 


_  ZB1  "  Zw 
XA  s 

w 


_  ZB3  ~  Zw 

XB  s 

w 


XC  XB 
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With  the  symbols  as  defined  in  Case  3,  and  h  =  -  z ^ 


AT  -  2  h  (xc  -  XA) 


XT  ' 


XA  +  XB  +  XC 


ZT 


2ZB1  +  ZB3 


AR  =  h<xBl  -  XC) 


-  _  XB1  +  XC 

XR  2 


ZB1  +  ZB3 


Then  the  submerged  area  Is 


the  force  of  buoyancy  is 


ZR 


As  =  AT  +  AR  ’ 


FB  =  P  AsWB 


FB 


'FB 


xtat  +  xrar 

A 


ZTAT  +  ZRAR 


A 


Case  5.  Both  rear  corners  (j  =  2,  j  =  4)  only  are  submerged 
(see  Figure  16) 


(34) 
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4(xB4’  zb4^ 


2^B2’  zB2^ 

FIGURE  16.  BUOYANCY  CALCULATION,  CASE  5 

TT 

Again,  since  9  is  never  greater  than  —  ,  C  can  be  assigned  without 
ambiguity  and  the  coordinates  of  A,  B  and  C  are  (x^,  z^) ,  (XB,  zBif) ,  and 
(x^,  z^),  respectively,  and: 


XA  = 


ZB2  "  Zw 


ZB4  "  Zw 


XB  S 


XC  XB 


With  the  symbols  defined  in  Cases  3  and  4, 


AT  2  h^xA  XC} 


XT  = 


XA  +  xB  +  xc 


2ZB2  +  ZB4 


AR  =  h<XC  "  xB2) 


XC  +  XB2 


XR  2 
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ZB2  +  ZB4 


Then 


A  =  A,  +  An 
s  T  R 


FB  =  P  As  wB 


x  jAy  +  x„A 


(35) 


R  R 


FB 


ztat  +  zRA( 


ZFB  A 


Case  6-  All  but  upper  front  corner  are  submerged  (see  Figure  17)* 

If  101  <  .01,  the  entire  box  is  assumed  submerged  and  the 
calculations  of  case  8  are  performed- 

3(*B3’  ZB3^ 


FIGURE  17-  BUOYANCY  CALCULATION,  CASE  6 

If  101  s  .01,  the  coordinates  of  A  and  B  are  (xn_,z.)  and  (xD,z__), 

B;  A  D  d3 


respectively,  where  z„  =  xD.s  +  z  and  xD  = 

r  '  AB3ww  B  s 

w 


The  symbols  are  as  defined  before,  and 


AT  2  ^ZA  "  ZB3^  ^XB3  "  V 
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2xB3  +  XB 


Zj  - 


2zB3  +  ZA 


Then  A  =  Ih 
s 


Aj  and 


XFB  ~ 


'FB 


FB  =  P  As  WB 


2  (XB1  +  xB2^h 


j(zBl  +  2B3)Ah 
A 


Vt 


"  ZT^T 


(36) 


Case  7*  All  but  upper  rear  corner  are  submerged  (see  Figure  18). 

If  1 9 1  <  .01,  it  is  assumed  that  the  entire  box  is  under 
water  and  the  results  are  as  in  Case  8,  below. 

4^XB4’  zb4^ 


If  101  ^  .01,  the  coordinates  of  A  and  B  are  (x^,  z^^)  and  (Xg^>  Zg) > 
respectively,  where  x. 


zD,  -  z 
B4  w 


A 


w 


and  zD  =  xD, s  +  z  • 
B  B4  w  w 


The  symbols  are  as  defined  before,  and 
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Then  A  =  -6h  -  A_  and 
s  T 


AT  2  ^ZB  "  ZB4^xA  "  xb4^ 


2xB4  +  XA 


-  _  2zB4  +  ZB 

ZT  3 


_  2  (XB 1  +  XB2^h  ~  xTAT 
XFB  “  A 

s 

2^ZB1  +  ZB3^h  ”  ZTAT 
ZFB  ”  A 

s 


(37) 


Case  8.  Entire  box  submerged. 

In  this  case,  with  the  symbols  defined  as  before, 

FB  =  P^h  WB 


FB  2^XB1  +  XB2^ 

(38) 

FB  =  2^Bl  +  ZB3^ 

Both  the  determination  of  which  case  exists,  and  the  calculations, 
are  coded  into  a  subroutine  called  BBOX,  described  in  Section  IV. 
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Equations  of  Motions 

The  equations  of  motion  used  here  are  a  restriction  of  those  derived 
by  Etkin^  for  six  degrees-of-f reedom,  body  centered  coordinates.  The 
restrictions  are  those  that  allow  only  vertical  plane  motion.  This  means 
that  surge,  roll  and  yaw  motions  are  assumed  constantly  zero. 

Let  n.  be  the  number  of  axles,  and  n„  be  the  number  of  body  boxes. 

A  D 

To  simplify  notation  in  this  section,  the  use  of  the  subscript  i  after  E 
implies  summation  over  all  axles;  i.e.,  i  =  l,,..,n^  ,  and  the  use  of  the 
subscript  j  implies  summation  over  all  body  boxes,  i.e.,  j  =  l,...,ng. 

The  surge  degree-of-f reedom  of  the  sprung  mass  is  governed  by  the 
equation 


M  (u-h*Q)  =  2  [E  F  .-sin  9  E  F,  .]  -  [W+E  FD.]sin  0  -  E  M  .a  .  (39) 
s'  '  sxi  bi  Bj  ui  xi 

Here  M  .a  .  is  the  force  exerted  on  the  body  due  to  the  acceleration  of 

thUI  X‘  12 

the  i  axle  in  the  x-direction  and  is  given  in  this  case  (see  Jurkat 

for  derivation)  by 

a^.  =  u  +  wQ  +  2Qz.  -  x.Q^  +  z.Q,  (40) 

Then,  substituting  a  .  into  Equation  (39)  and  solving  for  the  terms  in- 

X  I 

volving  derivatives 

Mu+(EM  .  z  .)Q.+2QE  M  .z.  =  -wQH  +2(SF  . -s  i  n  0E  F,  .  ]  -  [W+FD.]sin9 

'utr  uii  s  sxi  bi  Bj 

-  QEMu.(w-x.Q)  (41) 

Let  F  denote  the  right  hand  side  of  this  equation,  that  is 

X 

F  =  -wQM  +2  [EF  .-sin0  EF,  .]  -  [W+EF  .]s  i n 0-QEM  .(w-x.Q)  (42) 

The  heave  degree  of  freedom  is  governed  by  a  simpler  equation  since 
all  z-axis  motion  of  the  unsprung  masses  results  in  spring  and  damper 
forces,  not  inertial  forces  directly  on  the  sprung  mass.  Therefore,  in 
heave, 
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Ms(w  -  uQ)  =  -  ED.  +  (Ws  +  2  Fgj)  cos  0  (43) 

Again,  solving  for  the  terms  containing  derivatives,  and  renaming  the 
right-hand  side  Fz 

Mgw  =  MsuQ  -  ED.  +  (Wg  +  EFgj)cos9  =  F^  (44) 

The  pitch  equation  requires  organization  to  maintain  the  reader's 
orientation.  The  "applied"  moments  are  derived  from  the  following: 

o  suspension  spring  and  damper  forces 
o  wheel-soil  forces  in  the  x-direction 
o  weight  of  the  wheels  in  the  x-direction 
o  buoyancy  of  the  wheels  in  x-direction 
o  buoyancy  of  the  body  boxes  in  all  directions 
o  driving  torque  on  the  wheels 
o  tangential  soil  drag  on  the  wheels 
o  tangential  hydrodynamic  drag  on  the  wheels 
o  inertial  reactions 

The  pitch  equation  of  motion,  with  terms  written  in  the  above  order, 
then  becomes: 


I  Q  =  ED.x.  +  2ZM  .  -  EW  .z.  sin  9 
ys  ii  sxi  ui  i 


-  2EFb.Z.  sin  0  -  £FBj*F8j  cos0  '  FBj2FB  sin0 

-  T  +  E  (T.-D  .)  -  2  M  .z.a  . 

w  '  t  i  wr  uiixi 


After  substituting  a  .  from  Equation  (40)  into  Equation  (45)  and  solving 

X  I 

for  those  terms  containing  the  derivatives  (note  the  moments  of  inertia 

2 

of  the  unsprung  masses--  2Mu.z.  $--  add  to  IQ  to  produce  the  total  mass 
term  I  Q)  we  obtain: 

y 
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(EM  z  )u  +  I  Q  -  E(2QM  .z.)z.  =  SD.x.  +  2SM  . 
u  i  i  y  u  i  i  i  ii  sx i 

-  £W  .z.  sin  9  -  22  F,.z.  sin  9 
u  i  i  b  i  i 

'  EFBj  xFBj  005  6  '  EFBj  zFBj  sin  0  +  T„ 

+  E(TTi  '  Dwi)  -  d  EMu.2.(w  -  x.Q)  (46) 

Let  the  right  hand  side  of  Equation  (46)  be  denoted  by  Gy. 

The  equation  governing  the  rotational  speed  of  the  wheels  is 

I  q  =  T  +  2E  (T  .  -  D  .)  (47) 

The  equation  governing  the  position  of  each  wheel  center  relative 
to  the  vehicle  is 

M  .z.  =  W  .  cos  9  +  D.  +  2F  .  +  2FL .  cos  9  (48) 

uiiui  i  szi  bi 

Let  the  right  hand  side  of  Equation  (48)  be  denoted  by  Fz . . 

Required  for  the  simulation  are  the  position  and  velocity  of  the 
sprung-mass  center  of  gravity  in  bank  coordinates.  This  is  found  by 
calculating  the  velocity  of  the  sprung-mass  CG  in  bank  coordinates  from 

u  1  =  u  cos  9  +  w  s  in  9  (49) 

w1  =  -u  sin  9  +  w  cos  9  (50) 

and  integrating. 

The  equations  for  q  ,  u‘,  and  w‘  are  uncoupled  from  the  others  and 
may  be  treated  alone.  The  remaining  equations  are  coupled  and  cannot  be 
integrated  until  solved  as  a  system  for  the  derivative  terms.  To  this 
end,  we  write  them  as  a  linear  system  as  follows: 
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(51) 


Here  M 

u 

letting 


EH  .2.  and  the  convention  of  writing  the  wheel 

u  i  i  3 


equation  by 


2. 

I 


for 


1 . n 


A 


(52) 


was  followed. 


The  left-most  matrix  of  Equation  (51)  may  be  partitioned  along  the 
dashed  lines.  Therefore,  if  we  define  A,  B,  C  and  D  as  follows: 


A  = 


B  = 


I 


M 

0 

M 

L 

2QM 


0 

M 


u  1 


M 


M 


2QM 


u2 


..  2QM. 


un 


A 


2QM  ,2,  2QMu02  ...  2QM  z 
ul  I  22  unft  n^i 


(53) 
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C  =  (n.X3)  zero  matrix  =  0 

A  nAx3 

D  =  n.Xn.  identity  matrix  =  I 

A  A  1  nAxnA 

Then  the  system  of  Equations  (51)  may  be  written  as 

\ 


w 

Q 


z . 

I 


(53 


'A 


where 


J  = 


The  System  (51)  may  then  be  solved  by  finding  the  inverse  of  J. 

E  F 


J-1- 


H 


then 


E  =  A"1 
F  =  -  A-1B 


G  =  0 


V3 


H  =  I 


nAXnA 


This  means  that  only  the  3X3  matrix  A  needs  to  be  inverted.  Its 
is  given  by 


cont 1 d) 


(54) 


(55) 
If 

(56) 


(57) 


inverse 
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E  = 

M 


0  -M 


0  M/M 


-M  0 
u 


M  / 


(58) 


where  M  =  Ml  -  M^  .  Then 

y  u 


F  =  - 


2Q 


/ I  M  ,-M  M  .2. 
y  ul  u  ul  1 


I  M  -MM  z 
y  unA  u  unA  nA 


MM  ,z.-M  ,M 
ul  I  ul  u 


1 


MM  z  -M  M  , 
unA  nA  unA  u  j 


(59) 


Multiplying  the  matrix  J  thus  defined  by  the  vector  on  the  right  hand 
side  of  the  system  yields  the  following  set  of  first  order  equations, 
separated  by  degrees  of  freedom: 


Mu  =  |  F  -MG  -  2Q2  [M  .  (I  -M  z.)z  .] 
y  x  u  y  u  i  y  u  i  n^+ 1 

M  w  =  F 
s  z 

• 

MQ  --M  F+M  G  -  2QE  [M  .  (Mz.-M  )z  , 
u  x  y  u i  i  u  nA+i 


z.  =  z  . 
1  nA+ 


MuiZnA+l  Fzi 


I  D . 


(60) 

(61) 

(62) 

(63) 

(64) 


These  equations  may  now  be  integrated  separately. 
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IV.  SIMULATION 


The  simulation  described  by  the  previous  equations  is  mechanized 

in  a  FORTRAN  IV  program  written  for  the  Digital  Equipment  Corporation 

PDP-10  computer  located  in  the  Stevens  Institute  of  Technology  Computer 

Center.  The  program,  called  EGRES3-F4  is  run  under  the  time-sharing 

monitor  and  loads  into  13K  of  core-  Required  for  the  simulation,  in 

addition  to  various  subroutines  built  into  the  FORTRAN  compiler,  is  the 

1 2 

subroutine  HPCG,  part  of  the  Scientific  Subroutine  Package  (SSP)  ,  as¬ 
sembled  by  IBM  for  their  360  FORTRAN  compilers.  This  subroutine  controls 
the  simulation  and  performs  the  numerical  integration  of  the  equations  of 
motion  by  a  variable,  multi-step  fifth-order  predictor-corrector. 

Initially,  the  simulation  is  started  by  the  subroutine  INPUT  that 
first  reads  the  vehicle  and  bank  parameters,  which  are  stored  in  COMMON 
and  then  sets  the  initial  values  of  the  state  variables  and  stores  them 
in  an  array  called  Y.  The  derivatives  of  the  state  variables  are  then 
calculated  and  stored  in  an  array  called  DY  by  a  subroutine  called  UPDATE 
and  various  other  subroutines  called  by  UPDATE.  The  program  HPCG  then 
generates  "new"  values  of  the  state  variables  for  Y  and  calls  on  UPDATE 
to  compute  new  values  for  DY  as  a  result  of  tne  new  Y  values.  This  cycle 
keeps  repeating  until  the  simulation  is  terminated.  At  appropriate  times 
HPCG  calls  subroutine  OUTPUT  which  writes  the  values  of  the  state  variables 
in  Y  and  many  others  of  interest  onto  a  disk  file  for  future  printing. 

The  subroutine  hierarchy  is  given  in  Figure  19-  The  variables  shown  in 
Figure  19  are  only  the  significant  ones  passed  directly  from  subroutine  to 
subrouti ne • 
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Variables  in  common 


FIGURE  19-  ORGANIZATION  OF  EGRES3 -F4 
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The  storage  scheme  in  Y  is  as  follows: 

Y(l)  =  u  =  the  velocity  of  the  sprung  mass  CG  in  the  vehicle 

x-direction  (positive  is  forward) 

Y(2)  =  w  =  the  velocity  of  the  sprung  mass  CG  in  the  vehicle 

z-direction  (positive  is  down) 

Y(3)  -  0.  =  the  vehicle  pitch  velocity  (positive  is  nose  up) 

Y(4)  =  qw  =  the  rotational  velocity  of  the  wheels  (positive  is 

clockwise-rotation,  compatible  with  forward  motion) 

Y(5)  =  x^  =  the  x 1 -coord  inate  of  the  sprung  mass  CG  (positive  is 

to  right  of  the  water's  edge) 

Y(6)  =  z'  =  the  z '-coordinate  of  the  sprung  mass  CG  (positive  is 

down) 

Y (7)  =  0  =  the  pitch  angle  (positve  is  nose  up) 

Y(7+i)  =  zj  =  the  z-distance  from  vehicle  CG  of  the  center  of  wheel 

i  (positive  is  down) 

Y(7+n.+i)  =  z.  =  the  z-velocity  of  the  center  of  wheel  i  (positive 
is  down) 

Here  i  =  1 . n^.  Although  COMMON  contains  room  for  =  8,  Y,  DY  and 

AUX  are  currently  only  DIMENSIONED  large  enough  for  two  axles.  Thus  its 
size  must  be  changed  for  vehicles  with  more  than  two  axles. 

EGRES3.F4  requires  the  availability,  on  device  21,  of  a  file  con¬ 
taining  the  vehicle  parameters  and  a  file  containing  the  bank  soil  values 
and  profile.  These  files  are  to  be  in  the  format  described  in  Reference  3. 

The  following  describes  each  of  the  program  routines  in  more  detail. 
It  is  suggested  that  the  reader  have  before  him  the  listing  of  the  program 
and  follow  it  while  reading.  This  listing  may  be  found  in  Reference  5. 
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MAIN  Program 

The  MAIN  program  of  EGRES3.F4  serves  only  to  establish  storage,  to 
set  various  initializing  values,  to  call  INPUT,  and  then  to  turn  over 
control  of  the  simulation  to  HPCG.  Control  is  not  returned  to  MAIN  until 
the  simulation  is  to  be  terminated.  The  variable  IHLF  indicates  whether 
the  simulation  has  been  terminated  by  HPCG  due  to  numerical  difficulties 
(IHLF  =  11),  in  which  case  an  error  message  is  printed,  or  normally 
(IHLF  =  any  other  value). 

The  value  of  the  acceleration  due  to  gravity,  kept  in  GRAV,  sets 
the  units  of  the  simulation.  Currently,  GRAV  =  386.4,  which  implies  cal¬ 
culations  in  pounds,  inches  and  seconds. 

After  initializing  various  constants,  the  operator  is  asked  for  the 
name  of  the  output  file,  a  run  number  which  may  be  used  to  serialize  the 
various  simulations  made,  and  a  code  which  specifies  whether  and  what 
intermediate  output  is  to  be  printed.  (See  Reference  5  for  output  level 
code  specifications.)  MAIN  then  opens  the  output  file  and  passes  control 
to  INPUT. 

Upon  return  from  INPUT,  the  vehicle  parameters,  initial  positions 
and  velocities,  and  sufficient  bank  profile  points  to  span  twice  the 
vehicle  wheelbase  at  the  vehicles'  initial  position  are  in  storage.  MAIN 
then  transfers  the  initial  state  variables  into  the  array  Y  and  fills  DY 
with  the  error  bounds  used  by  HPCG  in  its  step  adjustment  procedures. 
Currently,  equal  bounds  in  all  degrees  of  freedom  are  used  if  a  two  axle 
vehicle  is  being  simulated.  Under  current  coding,  the  error  bounds  on  the 
wheel  positions  and  velocities  for  vehicles  with  more  than  two  axles  would 
be  more  stringent  than  the  body  degrees-of-f reedom.  This  may  not  be  the 
most  advantageous  weighing,  but  no  further  study  of  this  has  been  made. 

The  exact  way  of  using  these  error  bounds  is  given  in  the  description  of 

n 

HPCG  in  the  SSP  manual. 

After  initially  filling  Y  and  DY,  MAIN  sets  the  dimension  of  the 
equations  of  motion,  writes  the  first  output  page  title  and  headings  and 
transfers  control  to  HPCG. 
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Upon  return  from  HPCG  the  simulation  is  terminated. 

Subroutine  INPUT 

This  subroutine  reads  the  vehicle  parameters,  initial  values  and 
bank  profile.  INPUT  first  requests,  from  the  operator,  the  vehicle  data 
file  name.  INPUT  then  opens  that  file,  reads  it  and  echos  each  vehicle 
parameter  onto  the  output  file  opened  by  MAIN.  In  the  process,  component 
masses  are  calculated  from  their  weights. 

After  closing  the  vehicle  file,  the  name  of  the  bank  profile  file 
is  requested  from  the  operator.  INPUT  then  opens  that  file,  reads  it, 
echos  the  bank  title  and  soil  parameters  on  the  output  file,  and  then 
reads  XMIN(the  x-coordinate  of  the  first  bank  profile  point  given)5  DELTX 
(the  distance,  in  inches,  between  successive  bank  profile  points)^ and 
XMAX  (the  x'-coordinate  of  the  last  bank  profile  point  given). 

The  span  of  the  vehicle  is  next  calculated  by  adding  the  radius  of 
the  first  and  last  wheel  to  the  wheelbase  and  stored  in  XLTRTB.  This  is 
the  longest  extent  of  the  vehicle  which  can  touch  the  bank.  No  vehicle 
body-bank  interferences  are  considered;  these  should  be  added  for  future 
improvements . 

Since  the  numerical  integration,  to  be  accurate  and  stable,  "moves" 
the  vehicle  through  small  increments,  only  one  vehicle  span  need  be  in 
core  at  any  one  time,  as  long  as  the  bank  profile  directly  under  the 
vehicle  is  there.  In  order  not  to  access  to  disk  storage  too  frequently 
and  yet  not  to  use  too  much  of  core,  two  vehicle  spans  are  initially  read 
after  the  span  has  been  calculated.  XTMIN  contains  the  x'-coordinate  of 
the  leftmost  bank  profile  point  in  core,  XTMAX  the  x'-coordinate  of  the 
rightmost. 

The  array  BNK  contains  the  z'-coordinates  ( i .  e . ,  the  elevations)  of 
the  bank;  the  first  column,  the  original  bank  profile;  the  second  column, 
the  one  after  possible  modification  by  vehicle  travel.  Initially,  both 

*For  visualization,  the  vehicle  is  assumed  to  move  from  left  to  right  with 
the  bank  extending  from  the  lower  left  to  the  upper  right. 
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columns  are  filled  with  the  original  bank  profile  by  INPUT.  After  filling 
BNK  with  the  initial  two  vehicle  spans  of  the  bank,  INPUT  requests  the 
initial  values  from  the  operator.  Here  angles  are  converted  from  degrees 
(entered  by  the  operator)  to  radians  (used  by  the  simulation). 

Since  the  initial  position  of  the  vehicle  need  not  be  over  the  first 
two  vehicle  spans  of  the  bank,  the  x-coordinate  of  the  furthest  possible 
forward  part  of  the  vehicle  is  checked  against  XTMAX,  the  rightmost  bank 
profile  point  in  core.  If  the  bank  in  core  does  not  extend  over  two 
vehicle  spans  at  the  initial  vehicle  position,  more  bank  profile  data  is 
read  (10  points  at  a  time)  until  it  does.  If  the  end  of  the  BNK  storage 
array  is  reached  before  this  happens,  the  entire  bank  profile  is  shifted 
leftward  in  BNK  and  more  bank  read. 

After  BNK  contains  the  bank  profile  under  the  vehicle  at  its  initial 
position,  and  XTMAX  and  XTMIN  are  adjusted  appropriately,  the  number  of 
bank  profile  points,  KDTAU,  and  the  radius  of  the  largest  wheel,  DTAU,  is 
calculated.  Bank-wheel  interference  need  only  be  checked  for  a  distance 
DTAU  on  both  sides  of  a  wheel  center  x'-position. 

INPUT  then  returns  to  MAIN. 

Subroutine  UPDATE 

Subroutine  UPDATE  calculates  the  derivatives  of  the  state  variables 
whenever  called  upon  to  do  so  by  HPCG.  The  inputs  to  UPDATE  are  the  array 
Y  and  all  the  vehicle  and  bank  parameters  stored  in  COMMON.  The  output 
is  the  array  DY  containing: 


DY(1) 

• 

=  u 

DY(2) 

=  w 

DY  (3 ) 

• 

=  Q 

DY(4) 

DY(5) 

ii 

x» 

o  - 

CD 

II 

DY(6) 

=  *CG  = 

DY(7) 

=  e  =  Q 
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DY  (7+  i )  =  z  =  z  +i 

‘  nA 

DY(7+nA+l)  =  2. 

These  are  calculated  from  Equations  47  through  64.  This  is  a  rather  long 
and  complicated  procedure  encompassing  all  the  calculations  described  in 
the  section  entitled  EQUATIONS. 

Since  the  subroutine  is  called  at  least  twice  per  integration  step, 
some  concern  must  be  made  for  execution  time.  A  long  and  often  unnecessary 
step  in  FORTRAN  coding  is  the  repeated  use  of  subscripts.  In  UPDATE, 
whenever  a  subscripted  variable  must  be  used  more  than  twice,  it  has  been 
reassigned  to  a  non-subscr i pted  one.  In  fact,  the  first  action  of  UPDATE 
is  to  transfer  the  state  variables  from  the  first  seven  locations  of  Y, 
which  is  subscripted,  to  the  location  called  U,  W,  Q,  QW,  XCGP,  ZCGP,  and 
THETA.  The  parts  of  Y  describing  the  axle  motions  are  moved  to  subscripted 
variables  called  Z  and  ZD,  where  the  calculations  of  each  subscript  involves 
two  or  three  less  operations  than  the  subscripts  of  Y.  Later,  even  the 
values  in  Z  and  ZD  are  moved  to  unsubscri pted  variables  as  necessary. 

After  the  above  reassignment,  the  wheel  center  positions  in  the 
bank  coordinated  system  are  calculated  and  stored  in  XP(l)  and  Zp(|). 

The  corresponding  velocities  are  stored  in  UP ( I )  and  WP(|). 

UPDATE  then  checks  if  the  bank  profile  in  core  lies  under  the  current 
position  of  the  vehicle.  At  first  it  will,  since  INPUT  adjusted  it  so, 
but  after  some  time  the  vehicle  will  have  moved  rightward,  requiring 
additional  bank  profile  to  be  read  in  core.  The  procedure  here  is  similar 
to  that  in  INPUT.  In  case  the  vehicle  should  have  moved  left  beyond  the 
bank  profile  in  core,  an  error  message  is  printed  and  the  simulation  is 
terminated . 

The  next  major  task  of  UPDATE  is  the  calculation  of  the  entry  and 
exit  angles  for  each  wheel  with  respect  to  the  original  bank  profile  and 
the  current,  or  modified,  bank  profile.  A  test  for  wheel-bank  contact  is 
made  by  calculating  the  wheel  surface  position  directly  above  each  bank 
profile  point  for  the  distance  DTAU  on  each  side  of  the  wheel  center. 
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The  test  is  made  from  left  to  right  until  a  wheel  surface  point  below 
the  bank  profile  is  found.  If  none  is  found  until  a  distance  DTAU  on  the 
right  of  the  wheel  center  is  found,  no  wheel-bank  contact  is  made  and 
all  soil  forces  and  wheel-soil  variables  are  set  to  zero. 

The  first  wheel  surface  point  found  below  the  bank  elevation,  testing 
from  left  to  right,  gives  the  wheel  exit  angle  with  respect  to  the  parti¬ 
cular  bank.  If  the  original  bank  profile  is  used,  the  angle  so  found  is 
Y^o  °f  Figures  5  and  6  which  is  stored  in  GAM2(l,  1).  If  the  current 
bank  profile  is  used,  the  angle  found  is  and  is  stored  in  GAM2(I,  2). 

The  search  for  interference  is  then  started  from  the  right  going  to  the 
left.  The  first  wheel  surface  point  below  the  bank  elevation  thus  found 
yields  thewheel  entry  angle  Y^,  stored  in  GAM  1(1,  1),  or  Yj  stored  in 
GAM  1(1,  2).  In  all  these  cases,  wheel-bank  contact  is  assumed  midway  . 
between  the  bank  profile  point  below  the  wheel  surface  and  the  next  one 
above  the  wheel  surface.  Coding  is  included,  but  not  used,  which  allows 
more  precise  interpolation  if  required. 

The  location  of  the  wheel-bank  entrance  and  exit  points  are  then 
calculated  in  the  bank  coordinates,  and  stored  in  X1G,  Z1G  and  X2G,  Z2G. 
From  these,  the  local  ground  slope,  stored  in  ST,  is  calculated  and  then 
the  velocity  of  the  wheel  center  parallel  to  this  local  ground  slope  is 
calculated  and  stored  in  UPG.  Slip  is  then  calculated  by  Equation  7. 

The  entrance  and  exit  angles  with  respect  to  the  current  bank,  the 
entrance  angle  with  respect  to  the  original  bank,  the  ground  angle,  vehicle 
pitch  angle  and  the  state  variables  Y  are  now  passed  to  the  subroutine 
SOILF  and  the  values  of  T_  of  Equation  12,  F  of  Equation  17,  F  of 
Equation  18  and  M  of  Equation  19  are  returned  for  the  wheel  under  con- 

SX 

sideration. 

The  buoyant  forces  and  the  drag  on  the  rotating  wheel  are  calculated 
next.  In  the  program,  A  contains  the  area  of  the  entire  wheel  under  water, 
BB  the  part  of  the  rim  under  water.  Three  cases  are  distinguished: 

1)  the  entire  wheel  or  rim  out  of  water, 

2)  the  entire  wheel  or  rim  under  water,  and 

3)  part  of  wheel  or  rim  under  water. 
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The  equations  for  Cases  1)  and  2)  are  obvious.  In  Case  3)  the  area  sub¬ 
merged  is  calculated  using  Equation  20.  Wheel  buoyancy  is  calculated 
using  Equation  21  and  stored  in  FBW(l). 

To  calculate  rotational  drag  in  Case  1)  the  drag  is  set  to  zero; 
in  Case  2)  Equation  22  is  used  directly;  in  Case  3)  the  proportion  of  the 
circumference  under  water  is  applied  to  Equation  22.  The  rotational  drag 
of  the  I  wheel  is  stored  in  QWDRG(l). 

Engine  RPM  is  now  calculated  from  qw  for  each  gear  starting  at  low 
gear  and  going  to  MAXGR,  the  highest  gear  to  be  used  (requested  by  INPUT 
from  the  operator  at  initiation).  Negative  engine  RPM  results  in  an  error 
message  and  the  termination  of  the  simulation.  The  lowest  gear  yielding 
an  RPM  below  the  upper  shift  point  (in  UPSRPM)  is  used  as  the  operating 
gear.  If  no  gear  less  than  or  equal  to  MAXGR  yields  an  engine  RPM  below 
the  maximum  engine  RPM,  an  error  message  is  printed  and  the  simulation  is 
terminated.  The  engine  RPM  resulting  from  the  operating  gear  is  used  to 
determine  the  engine  torque  by  linear  interpolation.  Wheel  torque  is  then 
calculated  from  Equation  24  and  stored  in  WHLTRQ.  Since  it  is  assumed 
that  there  is  no  differential  rotation  among  the  wheels  and  all  wheels 
are  connected  to  the  one  engine  on  the  vehicle,  only  one  value  of  wheel 
torque  is  calculated  and  applied  to  all  wheels. 

The  suspension  forces  are  calculated  next.  The  location  of  each 

wheel  determines  which  of  the  Equations  25  is  to  be  used  and  the  spring 

force  is  calculated  and  stored  in  SPRGR.  Similarly,  a  test  is  made  to 

determine  which  of  the  Equations  26  is  to  be  used  for  the  damping  force; 

it  is  calculated  and  then  stored  in  DAMPF.  The  total  suspension  force 

t  h 

is  then  calculated  from  Equation  27  and  stored  in  D(l)  for  the  i  axle. 
The  values  used  to  describe  the  suspension  forces  are  the  sum  of  the 
forces  acting  on  both  sides  of  the  axle. 

The  buoyant  forces  on  the  body  are  calculated  by  subroutine  BBOX, 
described  below.  The  buoyant  force  on  box  J  is  stored  in  FBB(J) ;  the 
x-coordinate  of  the  center  of  buoyancy  of  box  J  is  stored  in  XCB(j)  ;  and 
the  z-coordinate  is  stored  in  ZCB(j). 
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All  the  values  required  for  the  equations  of  motion  have  now  been 
evaluated.  Appearing  in  the  equations  of  motion  are  16  summations  which 
are  evaluated  next.  These  are  stored  in  the  array  SM  according  to  the 
following  scheme: 

sm( l)  =  m  =  Em  .z. 

U  U  I  I 

SM(2)  =  £  M  . (I  -  M  z.)z. 

SM(3)  =EMu.(Mz.  -  Mu)z. 

SM(4)  =SFb. 

SH(5)=SFsxi 

SM(6)  =  E  Mu- (w  -  x.Q) 

SM(7)  =  Ed. 

SM(8)  =  E  D.  x. 

SM(9)  -  S«sxi 

SM(10)  =  E  (Tt.  -  D  .) 

SM(ll)  =  E  Ffa .  z. 

S M ( 12)  =  £W  .  z. 

u  i  i 

S M ( 13)  =£Mu.(w  -  x.Q)z. 

SM( 14)  =  E  FBj 
SH(15)  =SFBj  zfb. 

S"( 16)  -  S  Fej  xFB. 

Note:  In  storage,  z.  =  z^  +  .  ,  the  summation  over  i  goes  from  1  to  n^; 

A 

and  the  summation  over  j  goes  from  1  to  n^. 
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The  above  values,  and  others,  are  used  to  calculate  F  according  to 
Equation  42,  Fz  according  to  Equation  44  and  according  to  Equation  46. 
Then  the  derivatives  of  the  state  variables  are  calculated  according  to 
Equations  60,  61,  62,  47,  49  and  50.  The  derivatives  of  the  wheel 
position  and  velocity  are  calculated  from  Equations  63,  64  and  48. 

Subroutine  OUTPUT 

Subroutine  output,  called  by  HPCG  whenever  an  integration  step  is 
completed,  determines  if  sufficient  simulation  time  has  passed  to  print 
an  output  or  sufficient  distance  has  been  traveled  to  update  the  bank 
profile.  If  so,  it  takes  the  appropriate  actions. 

Since  there  is  a  variable  number  of  lines  in  the  output  block 
'(depending  on  the  number  of  axles,  n^  ,  and  the  number  of  body  boxes, 
n^)  OUTPUT  also  checks  if  enough  space  is  left  on  the  current  output 
page  for  another  output  block.  If  there  is,  it  prints  it;  if  not,  it 
starts  a  new  page  and  prints  the  run  number,  bank  title  and  output  block 
legend  before  the  next  output  block. 

Besides  reporting  on  the  progress  of  the  simulation  whenever 
necessary,  OUTPUT  also  modifies  the  bank  to  conform  to  wheel  travel. 

This  cannot  be  done  in  UPDATE,  since  the  equations  of  motion  are  often 
calculated  for  trial  values  of  the  integration  interval.  If,  as  a 
result  of  a  trial  value,  the  error  made  between  the  predictor  and 
corrector  (in  HPCG)  is  too  large  or  too  small,  the  integration  interval 
is  modified  and  the  equations  of  motion  recalculated.  If  the  bank  had 
been  modified  in  UPDATE,  it  would  have  to  be  restored  when  false  inte¬ 
gration  steps  are  used.  This  would  have  been  very  complicated,  and 
unnecessary. 

HPCG  continues  to  modify  time  increments  until  a  successful  inte¬ 
gration  step  has  been  achieved.  It  then  transfers  control  to  OUTPUT. 

Now  the  bank  profile  can  be  modified  to  conform  to  the  current  wheel 
positions.  This  is  done  by  determining  all  bank  profile  points  at 
which  bank-wheel  interference  occurs,  and  setting  the  bank  elevation 
at  the  elevation  of  the  wheel  edge.  This  is  done  only  with  the  current 
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bank  profile  stored  in  BNK(1,  2).  The  original  bank,  in  BNK ( 1 ,  1)  is 
left  unmodified. 

Subroutine  SOILF 

This  subroutine,  called  for  each  axle,  calculates  the  forces  on 
the  wheel  and  body  due  to  wheel -soil  interference.  This  coding  is  for 
the  rigid  wheel  model  described  in  Misklevitz.  The  wheel  is  specified 
by  the  first  variable  in  the  calling  sequence. 

First,  SOILF  determines  the  appropriate  soil  parameters  to  be  used 
for  the  current  position  of  the  wheel.  It  then  calculates  the  maximum 
sinkage  by  Equation  4,  and  maximum  normal  stress  by  Equation  3*  The 
coefficients  (Equation  2)  of  the  normal  stress  distribution  given  by 
Equation  I  are  then  calculated. 

The  integrals  in  Equations  8  and  9  are  integrated  analytically  by 
Misklevitz. ^  The  resulting  formulas  are  used  to  calculate  for  the  I ^ 
wheel  F  ,  stored  in  SIGH(l),  and  F  ,  stored  in  SIGN(l). 

Since  the  integrals  of  Equations  10  through  14  were  not  able  to  be 

14 

evaluated  analytically,  a  five  point  Newton-Coles  formula  is  used  to 
evaluate  them.  The  shear  stresses  at  five  angles  between  Yj  and  was 
evaluated  using  Equations  5  and  6.  Figure  20  shows  the  location  of  the 
five  angles  chosen  for  the  numerical  integration.  These  five  angles  were 
picked  after  considerations  given  by  Misklevitz. ^  The  five  point  Newton-Cotes 
formula  then  gave  F^,  stored  in  TAUN(l);  F^T,  stored  in  TAUH(l);  and  , 
stored  in  TTAU(  I )  . 

F  (Equation  17)  and  F  (Equation  18)  are  then  calculated, 
s  X  s  z 

Though  somewhat  more  complicated,  M  was  calculated  by  a  similar 

sx 

Newton-Cotes  formula- 
Subroutine  BBOX 

This  subroutine  calculates  the  buoyant  forces  and  center  of  buoy¬ 
ancy  of  the  vehicle  body  box  specified  by  the  first  variable  in  the 
calling  sequence . 
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After  transferring  the  vehicle  coordinates  of  the  box  corner  to 
unindexed  variables,  their  bank  coordinates  are  calculated  using 
Equation  29.  A  sequence  of  tests  is  then  made  to  determine  which  of 
the  cases  described  in  section  entitled  Body  Buoyant  Forces  applies. 

The  remainder  of  the  coding  of  this  subroutine  follows  the  equa¬ 
tions  of  that  section  so  closely  that  no  further  explanation  is  needed. 


FIGURE  20.  TYPICAL  SHEAR  STRESS  DISTRIBUTION 
SHOWING  POINTS  USED  IN  NUMERICAL 
APPROXIMATION  OF  ITS  INTEGRAL 


T  (y) 
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NOMENCLATURE 

Unprimed  variables  measured  in  vehicle  coordinate. 

Primed  variables  measured  in  earth-fixed  coordinate. 

Dots  above  variables  mean  differentiation  with  respect  to  time. 

Subscripts 

b  wheel  buoyancy 

B  refers  to  body  buoyancy 

c  wheel  compression  (or  jounce)  direction 

CG  center  of  gravity  of  sprung  mass 

G  parallel  to  local  ground  slope  or  refers  to  transmission 

i  wheel  index,  i  =  l,...n^  ;  gear  index  i  =  1,...,N^ 

d  damper 

p  spring 

r  wheel  rebound  direction 

R  rim,  wheel  rest  position,  engine  speed,  rectangle 

T  triangle 

s  soil,  submerged,  or  sprung  mass 

u  upper  suspension  bump  stop,  or  unsprung  mass 

x  in  vehicle  x-direction 

z  in  vehicle  z-direction 

Y  at  point  on  wheel  surface  located  by  a  line  from  wheel 

center  making  an  angle  y  with  respect  to  horizontal 

j  body  box  index,  j  =  l,...,ng;  or  box  corner  index,  j  =  1,2,3,^ 

i  lower  suspension  bump  stop 

cu  wheel 

a  arising  from  soil -wheel  normal  force 

t  arising  from  soil -wheel  shear  force 
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Variables 
a  ,b 
A  u 

3  D  C  •  •  • 

A,B,C 


b 

c,cp 

d  ,d 
r’  c 

D 

Dw 

ERi 


abc. • • 

F,  ,F. 
dr  dc 


G. 


I  >1 

y  ys 

j 

k 

k  ,  k  ,n 
c’  cp 

k  ,k 

u  l 


M 


H 

abc*  •  • 
M 

u 


n 

n 

n 


B 

A 

G 


Nomencl ature 
(cont'd) 

parameters  of  wheel  rotational  drag  equations 
area  designated  by  subscripts  a,b,c,.*. 

points  used  in  vehicle  body  buoyancy  locations  or 
coefficients  of  normal  stress  equations 

wheel  width 

Coulomb  soil  shear  strength  parameter 
wheel  damper  (shock  absorber)  constant 

total  suspension  force 
wheel  rotational  drag 
,  .  .  th 

engine  speed  in  i  gear 

force  designated  by  subscripts  a,b,c,**. 

maximum  force  derivable  from  wheel  dampers  (blow-off  force) 

final  drive  ratio 
.*  r  •  th 

gear  ratio  of  i  gear 
total  appl  ied  moment  about  CG 
height  of  box 

pitch  moment  of  inertia  of  total  vehicle  and  sprung  mass 
soil  deformation 

spring  constant  or  soil  deformation  parameter 
Bekker's  soil  strength  parameters 

spring  constant  of  upper  and  lower  bump  stops 

total  vehicle  mass 

moments  or  masses  designated  by  subscripts  a,b,c,*»* 

total  mass  moment  of  unsprung  masses 

number  of  body  boxes 

number  of  axles 

number  of  transmission  gears 
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Nomencl  ature 
(cont  'd) 


Var iab les 


P. 

i 


w 


r 

s 

S. 


S  5  S  t 

P  a 


W 


abc  •  •  • 


oo 

u 


w 

wr 


FB 


'FB 


z . ,  z . 

i  i 


bank  profile  location  points  (in  x'-z1  coordinates) 

vehicle  pitch  velocity  (positive  nose  rising) 

wheel  rotational  velocity  (positive  when  wheel  bottom 
moving  towards  rear  of  vehicle) 

wheel  radius 

wheel  rotational  slip 

bank  x 1 -coord  ina  tes  at  which  one  or  more  soil  parameters  change 

spring  force  due  to  spring  and  damper 

engine  speed  at  which  next  lower  gear  is  selected 

slope  of  waterline  with  respect  to  vehicle  coordinate 

torques  designated  by  subscripts  a,b,c,... 

engine  torque 

torque  at  wheels 

vehicle  surge  velocity  (positive  forward) 

velocity  of  wheel  center  parallel  to  local  ground  slope 

engine  speed  at  which  next  higher  gear  is  selected 

vehicle  heave  velocity  (positive  down) 
width  of  body  buoyancy  boxes 

surge  direction  moving  with  vehicle  (positive  forward) 

bank  coordinate  parallel  to  water  surface  (positive  bankward 
from  water  bank  intersection) 

x-coordinate  of  box  center  of  buoyancy 

heave  direction  moving  with  vehicle  (positive  down) 

bank  coordinate  perpendicular  to  water  surface 
(positive  down  for  water  surface) 

z-coordinate  of  box  center  of  buoyancy 
location  of  center  of  wheel  i 
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Nomenclature 

(cont'd) 


wheel  center  position  when  lower  bump  stop  is  first  contacted 
wheel  sinkage  of  center  of  contact  patch 

wheel  center  position  when  spring  is  extended  to  its  no-force 
distance 

wheel  center  position  when  upper  bump  stop  is  first  contacted 
z-intercept  of  waterline  in  vehicle  coordinates 

angle  defined  by  tt/2  +  y  -  ft 

efficiencies  of  transmission  and  final  drive 

angle  made  by  line  connecting  wheel  center  with  a  point 
in  the  contact  patch 

entry  angle  of  wheel  with  respect  to  modified  bank 

exit  angle  of  wheel  with  respect  to  modified  bank 

entry  angle  of  wheel  with  respect  to  original  bank 

exit  angle  of  wheel  with  respect  to  original  bank 

length  of  box 
dens i ty  of  water 
soil -wheel  normal  stress 
soil -wheel  shear  stress 

pitch  rotational  translation  moving  with  vehicle  (positive 
nose  up  measured  from  x-axis) 

local  bank  slope 

angle  with  respect  to  vertical  made  by  line  from  wheel  center 
to  wheel-water  intersection 
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